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Abstract 

We modify the time dependent Schrodinger-Newton equation by using a poten- 
tial for a solid sphere suggested by Jaaskelainen (Jaaskelainen 2012 Phys. Rev. 
A 86 052105). Compared to our recent paper (Giulini and GroBardt 2011 Class. 
Quantum Grav. 28 195026) where a single point-particle, i. e. a Coulomb po- 
tential, was considered this has been suggested to be a more realistic model for a 
molecule. Surprisingly, compared to our previous results, inhibitions of disper- 
sion of a Gaussian wave packet occur at even smaller masses for the solid sphere 
potential, given that the width of the wave packet is not exceeded by the radius 
of the sphere. 



1 Introduction 

In UJ we investigated the time dependent Schrodinger-Newton equation (henceforth 
abbreviated SN-equation): 

ihdt*(t,x) = (-— A-Gm 2 f Ifi*'^), 2 d 3 y] *(t,x). (1) 
V 2m J \\x-y\\ ) 

The potential term can be thought of as the convolution of the absolute value squared 
of the wave function with the Coulomb potential 3>(r) = —Gm 2 /r: 

-Gm 2 f ^|f ^ = (^|$p)(t,g), (2) 

J If- 2/11 

Our numerical and analytical consideration of the SN equation in 111 showed that 
inhibitions of dispersion due to the gravitational self-interaction occur at mass values 
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of about 6.5 x 10 9 u for a given wave packet width of 0.5 fim.. This result was contrary 
to a previous one by Salzman and Carlip (2l|3l, where inhibitions of dispersion were 
claimed to occur at mass values more than six orders of magnitude smaller, namely at 
m m 1600 u. Somewhat loosely, this behaviour of a wave packet initially shrinking in 
width was called a "collapse" by Salzman and Carlip an we adopt this nomenclature. 

As we pointed out in the Summary of @, more consideration has to be given to 
the question concerning the application of the SN equation to realistic systems, like 
molecules. In a recent paper Q Jaaskelainen suggested that the gravitational self- 
interaction for a molecule should be modelled by the potential of a solid sphere of 
radius R 

( Gm 2 ( 3 - 2 



$(r) = t R V 2 2&J> lir<R (3) 
I Hr>R 
which has to be put into equation Q instead of the Coulomb potential. 

Jaaskelainen points out that in our case, for a wave packet width of 0.5 fim, the 
results obtained for this potential deviate from the results for the Coulomb potential 
for masses beyond 10 10 u. 

We will not make any statement about the validity of this ansatz but consider the 
deviations that one obtains, if the Coulomb potential in our previous analysis in [1] is 
replaced by the solid sphere potential ([3]). 



2 Results 

We study the evolution of the SN equation making use of the same numerical methods 
used in HI, simply replacing the potential by the solid sphere potential Q. All results 
refer to an initial Gaussian wave packet 

*( r ,t = 0) = (vra 2 )" 3 / 4 exp (~^) (4) 

with a width of a = 0.5 /im. 

First we look at the evolution of wave packets if the density of the solid sphere 
corresponds to the density of either lithium (534 kg m~ 3 ) or osmium (22610 kg m" 3 ), 
cf. 0. To illustrate the results, the radius within which 90 % of the probability density 
are contained is plotted against time in figure [I] for several masses. In addition to the 
results for the density of lithium (Li) and osmium (Os) the results for the SN equation 
with Coulomb potential obtained in [1] are also plotted. 

For all masses considered the modified potential results in an even faster collapsing 
wave packet. The shrinking behaviour also sets in for smaller masses, namely for 
5 x 10 9 u the modified solution shrinks while the SN equation solution is still spreading 
for this mass. 

To get a better impression of this behaviour, let us take a look at the effect that the 
radius of the solid sphere potential has on the mass for which the collapsing behaviour 
of the wave packet sets in. We thus ask the following question: given a sphere of radius 
R, what mass m c do we have to distribute homogeneously "within" this sphereQfor a 
collapse of the wave packet to take place? 

1 But note that the gravitational interaction is always weighted with the probability density of the Gaus- 
sian wave packet whose width is always fixed to 0.5 fim. 
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Figure 1: This plot shows the radius rgo within which 90 % of the probability lie plotted 
against time for several masses for the SN equation with Coulomb potential (SN) and with 
a solid sphere potential for both the density of lithium (Li) and osmium (Os). 



In HI we found that for a 0.5 /mi wave packet a collapse is observed for masses of 
about 6.5 x 10 9 u and higher. In figure [2] this critical mass value is plotted depending 
on the radius of the solid sphere. 

For small radii of several nanometers this critical mass reduces to about two-thirds 
of the value obtained in [ 1 ] and it stays smaller than this value until the size of the 
solid sphere exceeds the width of the wave packet significantly. There seems to be an 
extremal value somewhere in between 20 and 100 nm. 



2.1 Why does R > have a positive effect on the collapse mass? 

In section 3.2 of [T] we derived a total energy for the SN wave packet: 

h 2 8Gm 2 



£ 
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-I(R/a) 



with the integral expression 



1(0) 
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(5) 



Only initial data that has a negative total energy can show inhibitions of dispersion (6l 
and therefore the zero of the total energy ([5]) provides an estimate value for the collapse 
mass. 
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Figure 2: The critical mass beyond which the wave packet reveals a shrinking behaviour 
is plotted against the radius of the solid sphere in the potential term. 



Let us see what changes if we modify the SN equation employing the solid sphere 
potential (|3]>. The total energy is unchanged except of the integral expression, which 
now takes the following form: 
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(6) 
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Figure 3: Plot of the integral expression I(R). 



As I(R) is strictly positive, inhibitions of dispersion can only occur for masses 



m > 



ir ft 2 a 



16 GI(R/a) 



1/3 



(7) 



Thus, the minimal collapse mass can be obtained, where I(R) reaches its maximum. 
The function is plotted in figure [3] The maximum value is taken for R rj 0.598 a. 
Although this analysis only yields a lower mass bound, this value fits well to the nu- 
merical results described previously. 



3 Conclusion 

If one does not take the simpler point of view we took in [ 1 ] and model large molecules 
by a point-particle SN equation, and instead assumes the modified equation given in 
to provide a more realistic description for such molecules, our results confirm that 
there are indeed deviations between both equations, as Jaaskelainen suggested. 

However, whereas one naively might expect that a smaller matter density leads 
to a weaker gravitational interaction and therefore to a behaviour that is closer to the 
free Schrodinger equation, our results show that the converse is true. As long as the 
radius of the solid sphere potential does not significantly exceed the extent of the wave 
packet, the effect of inhibitions of dispersion is even increased compared to the SN 
equation. The analysis of the total energy given in section [2T] provides an insight why 
this somewhat peculiar behaviour shows up. 

The question why there should be a self-gravitational interaction while there is no 
electromagnetic self-interaction term present in the Schrodinger equation remains an 
issue that should be addressed (cf. the discussion in O). This issue arises even more 
distinctly in the analysis given in [5] where the mutual gravitational interactions of the 
constituents of a molecule are assumed to be negligible compared to the electromag- 
netic ones while gravitational self-interaction has to be considered only because there 
is no electromagnetic ^//'-interaction. 

We conclude that taking into account the finite extent of the source in the way 
suggested in JH does not lead to an attenuation of the effect of inhibitions of dispersion 
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as discussed in [1]. The prospects for an experimental verification of self-gravitation of 
quantum systems in future molecular interferometry experiments remain unchanged. 
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